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(—1 Abstract 

^"The IceCube experiment at South Pole consists of two detector components - the IceTop air shower array on the surface and the 
Q neutrino telescope at depths from 1450 to 2450 meters below. Currently, 26 IceTop stations and 22 Inlce strings are deployed. With 
^— < the present size of the IceTop array, it is possible to measure cosmic rays with energies ranging from 0.5 PeV to 100 PeV. Coincident 
CZJ events between the IceTop and the Inlce detector provide useful cross-checks of the detector performance and furthermore make 
it possible to study the cosmic-ray composition. This paper gives an overview on the current status of IceTop. 
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1. Introduction 

The IceCube neutrino telescope [T] being built at South 
Pole, consists of two major detector components. The ac- 
tual neutrino telescope (Inlce) is buried at depths between 
1450 m and 2450 m in the antarctic ice and will con- 
sist of 4800 individual photon detectors (Digital Optical 
Modules), lined up on 80 strings in an hexagonal pattern. 
The IceTop air shower array is located on the ice surface 
above the Inlce detector. It consists of 80 stations, close to 
the position of the Inlce strings, with a spacing of 125 m. 
Presently 26 stations are operational and are arranged as 
shown in FigfTJ 

The IceTop air shower array is a multi-purpose detector 
with several scientific aspects and technical goals. It helps 
to understand the background of single muons in the In- 
lce detector and acts as a veto against muon bundles from 
cosmic ray induced extensive air showers (EAS). Further- 
more, it is used independently to detect high energy cosmic 
rays with energies ranging from approximately 500 TeV up 
to about 1 EeV (in its final stage). Coincident air shower 
events, measured by the IceTop and the Inlce detector, 
make it possible to study the cosmic-ray composition since 
the number of muons, which are able to reach the Inlce de- 
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tector, is also sensitive to the mass of the primary cosmic- 
ray particle. In addition, coincident events provide very 
useful cross-checks for the precision of the direction- and 
energy-reconstruction in IceTop and Inlce. The following 
sections summarize the current status of IceTop and give 
some examples for the performance of the detector. 

2. The IceTop Array 

Each of the 80 IceTop stations consists of two Ice- 
Cherenkov tanks with an inner diameter of 1.86 m and an 
ice thickness of 0.9 m as shown in Fig. [5] The tanks are 
buried under a thin layer of snow at a distance of about 
10 m apart. The inner surface of the tanks is covered with 
a diffusive coating to homogeneously reflect the induced 
Cherenkov light. Two DOMs, identical to those used in the 
Inlce detector, are frozen half submerged in the ice surface 
and measure the Cherenkov emissions in the tank. One 
DOM is operated at high gain (HG) and the other at low 
gain (LG) in order to extend the dynamic range of the tank. 

After deployment of a tank the freezing process of the 
ice is controlled by a freeze control unit (FCU) which circu- 
lates and degases the water in the tank while the temper- 
ature is slowly decreased. This way the water freezes from 
the top to the bottom of the tank and the number of bub- 
bles is significantly reduced [2]. The freezing process of a 
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Fig. 1. Event display with the present layout of the 26 station IceTop 
array. Each station consists of two tanks, represented by a half circle. 
The radius of the half circles scales with the signal in the tanks and 
the different colors correspond to the signal times. The reconstructed 
core and direction are illustrated by the star and the arrow. The 
dashed line traces the shower plane through the core on the surface 
of the array. 

tank usually takes six to eight weeks. After the freezing is 
completed the FCUs are recovered and refurbished to use 
them in the next deployment season. 

The two tanks that make up a station are linked together 
by local coincidence conditions between their DOMs in or- 
der to reduce the number of accidental coincidences. A lo- 
cal coincidence basically occurs if two DOMs in different 
tanks of a station measure a signal within 250 ns. If at least 
six IceTop DOMs report a local coincidence within 2 /is 
a Simple Multiplicity Trigger (SMT) is generated, causing 
the data acquisition to read out the air shower event. The 
SMT trigger rate depends on the size of the array and is 
currently of the order of 14 Hz. In addition the IceTop ar- 
ray is always read out when the Inlce detector has been 
triggered or vice versa. 

The IceTop tanks are calibrated conventionally using at- 
mospheric muons to determine the mean signal of a Ver- 
tical Equivalent Muon (VEM) in the tanks. This unit can 
be easily compared to air shower simulations. Muon tele- 
scope measurements and simulations show that the actual 
muon peak is at about 95 % of the peak position of the full 
recorded spectrum of all hits [3] . 
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Fig. 2. Sketch of an IceTop tank containing a high gain (HG) and a 
low gain (LG) DOM. The inner surface of the tanks is covered by a 
reflective coating. 

of the reconstruction and the uncertainties of the various 
shower observables. An example for the Arc-distribution of 
the core reconstruction for the 16 station array in 2006 is 
given in Fig. [3] where only contained events which fulfill the 
SMT trigger condition were taken into account. The dis- 
tribution in Fig. [3] is fitted by a Gaussian function. Since 
the spread <jax of this distribution emerges from the fluc- 
tuations of Ax = xa — xb the uncertainties oa,b of the 
individually reconstructed x-coordinates xa ot xb can be 
obtained from 
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Assuming that the actual value of x, using the complete ar- 
ray, can be approximated as the average (x) — (xa + xb)/2 
of the results of the two sub-arrays, the actual uncertainty 
of (x) becomes 
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which is half of the spread of the Ax-distribution shown 
in Fig. [3l The same procedure was applied to the other 
shower observables revealing a core reconstruction preci- 
sion of about 17 m, and an uncertainty in the zenith and 
azimuth angles of 1° and 2° for the 16 station array in 2006. 
The reconstruction quality will increase further with the 
size of the IceTop array. Since the results of this analysis 
may depend on the reconstruction algorithm, the obtained 
uncertainties shouldn't be interpreted as detector resolu- 
tions but are a valuable cross-check for the stability of the 
reconstructions. 



3. Reconstruction Accuracy &; Stability 

Due to the special layout of the IceTop array with two 
tanks (A and B) per station, it is possible to subdivide the 
array into two identical sub-arrays (A and B) with just one 
tank per station. This way an air shower event can be inde- 
pendently reconstructed by the two sub-arrays. The com- 
parison of the results gives an estimate for the stability 



4. Single Station Coincidences 

A precise time synchronization among all DOMs and a 
good knowledge about their position is a critical require- 
ment for any physics analysis. Calibration with flashers and 
survey by hole logging during deployment shows that the 
time synchronization of an Inlce string is at the level of 3 ns 
while the depth of the individual DOMs are known with an 
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Fig. 3. Ax-distribution of the individually reconstructed x-coordi- 
nates of the shower cores in the two sub-arrays. See text for details. 

accuracy of about 50 cm pQ . The timing between the Ice- 
Top and the Inlce detector components can be checked by 
using vertical muons which were tagged by triggering only 
a single station in the IceTop detector. To ensure that the 
single station events are not caused by tails of big air show- 
ers outside the array, only the inner stations of the IceTop 
array are used together with the Inlce strings directly be- 
low them. With the 16 IceTop stations and 9 Inlce strings 
in 2006, only stations 39 and 49 fulfill this requirement. For 
these two strings the muon speed has been individually cal- 
culated for each Inlce DOM relative to the time to of the 
IceTop station at the surface according to u, = dk/ (i, — to) 
where di is the distance between the station and the i th 
Inlce DOM [4]. In first order the arrival times of the pho- 
tons at the DOMs are exponentially distributed due to the 
scattering in the ice. Since there are a number of other in- 
dependent uncertainties in the system the exponential dis- 
tribution is smeared out according to a Gaussian function. 
Therefore the time distributions of individual DOMs have 
been fitted with a Gaussian convoluted exponential func- 
tion 
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where the time constant r corresponds to the scattering 
length in the ice and a is the effective time resolution of 
the Gaussian function. The expression "erfc" represents the 
complementary error function [5]. An example of such a 
fit is given by the graph included in Fig. |U The fitted ti 
values for each DOM correspond to the time origin of the 
pure exponential function and thus to the time offset of 
the DOM. The distribution of muon speeds, in units of the 
speed of light c, was determined using these time offsets 
and is shown in Fig|4] where the data of the two strings 
has been combined. Only the data of 110 out of available 
120 DOMs were used since the fit failed in the other cases 
due to insufficient statistics. As can be seen in FigfJ] the 
muon speed in the ice is compatible with the speed of light 
and the RMS of 0.0015 reflects the uncertainties in the 
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Fig. 4. The distribution of muon speed (v) relative to the speed of 
light (c) [4]. The cut-in entry shows the time delay on one in-ice 
DOM and the fit. See text for details. 

timing and the location of the DOMs and of the true muon 
position on the surface. If only the timing uncertainty is 
taken into account the RMS translates to an upper limit 
on the timing precision of about 12 ns in 2.5 km depths. 
Although the precision of this method is not as good as the 
standard survey and calibration techniques, it is a useful 
and independent cross-check that there is no significant 
deviation from expectation. 



5. Energy Spectrum &: Composition 

The full reconstruction of an IceTop event is an itera- 
tive process. First, the center of gravity (COG) of all tank 
signals is used as a first guess for the core position. In a 
second step the arrival direction of the air shower is re- 
constructed, using the time information of the signals and 
assuming the shower front to be a plane wave. These first 
guess values are taken as the seeds for a final log-likelihood 
minimization which fits a lateral distribution function to 
the tank signals calibrated in VEM. For IceTop a Double 
Logarithmic Parabola (DLP) was found to describe the lat- 
eral distribution the best [6]. The reconstructed signal at 
a distance of 100 m from the shower axis (Siqo) is used as 
an estimator for the primary cosmic-ray energy. The con- 
version between the S'ioo value and the energy is based on 
a nearly linear relation between the two quantities but also 
depends on the zenith angle of the shower [6] . For showers 
with zenith angles less than 30° the mean energy for S'ioo = 
20 VEM is approximately 10 PeV and S'ioo = 200 VEM 
corresponds to roughly 100 PeV [7]. The resulting raw en- 
ergy spectrum without any acceptance correction is shown 
in Fig. [5] At high energies, where the events are detected 
and reconstructed with high efficiency, the slope of the spec- 
trum agrees well with the spectral index 7 = 3.05 found by 
other experiments |6j which is represented by the solid line 
in Fig. O 

Coincident air shower events, detected by IceTop and In- 
lce, make it possible to determine the primary energy in 
three different ways, making IceCube a hybrid cosmic-ray 
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Fig. 5. Preliminary raw energy spectrum without acceptance correc- 
tion. In the high energy range the spectrum already agrees well with 
the expected flux from other experiments [6]. 

detector and providing a cross-check for the different re- 
construction methods and their systematic uncertainties. 
Besides the method expiained above, also the lateral dis- 
tribution of Cherenkov photons from the muon bundles in 
the Inlce detector can be used to estimate the primary en- 
ergy. This reconstruction is currently being developed and 
provides an independent approach to estimate the primary 
energy, however with a coarse energy resolution due to the 
primary-mass sensitivity of the muon number. Also the ge- 
ometry of the shower can be independently reconstructed 
by the Inlce detector. In last instance a coincident event 
can further be reconstructed by combining the IceTop and 
the Inlce data. This method leads to precise and consistent 
results in the whole IceCube detector. 

Since the number of muons in an air shower depends 
on the energy and the nature of the primary particle, the 
cosmic-ray composition can be studied by correlating the 
muon and the electron numbers as is suggested by the sim- 
ulation study shown in Fig. [6l The energy increases along 
the diagonal of this plot whereas the mass of the primary 
goes from light (proton) to heavy (iron), roughly perpen- 
dicular to the hypothetical energy axis. However, in real 
data these scatter-plots look much more fuzzy due to the 
limited detector resolution and the occurrence of all kinds 
of primary particles from proton to iron. The individual 
energy spectra for different kinds of elemental groups can 
be statistically disentangled with different techniques, for 
instance by a two dimensional unfolding approach as used 
in [8]. The muon reconstruction code for IceCube/IceTop is 
still under development. The cosmic-ray composition study 
is one of the most important scientific goals for IceTop. 

6. Summary 

The IceTop air shower array at South Pole grows each 
austral summer. After the coming deployment season, 50 % 
of the array will be completed. In parallel to the detector 
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Fig. 6. Simulation of the correlation between the muon and the 
electron numbers in an air shower [9]. 

construction the reconstruction tools are being developed. 
An overview of the detector performance and the existing 
or planned scientific reconstruction methods has been given 
in this article. The current reconstruction effort focuses on 
the development of a new muon bundle reconstruction for 
IceCube to improve the energy resolution for coincident air 
shower events and to study the cosmic-ray composition. 
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